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Abstract. Although the formulas for the light deflection due to quadrupole 
gravitational field of deflecting bodies are well known, the formulas are rather 
complicated, so that massive computations of quadrupole light deflection (e.g., in the 
framework of astrometric survey missions like Gaia) are time-consuming. Gonsidering 
an observer situated within a few million kilometers from the Earth (clearly the most 
practical case), we derive the simplest possible form of the relevant formulas still 
having numerical accuracy of 1 /zas. This form leads to simple upper estimates for 
the quadrupole light deflection in various cases allowing one to relate the magnitude of 
the actual quadrupole deflection with the corresponding monopole deflection due to the 
same body. These upper estimates can be used to decide if, for a given conflguration, 
the actual quadrupole deflection should be computed for a given accuracy goal. 
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1. Introduction 

In the nearest future astrometric observations will reach an accuracy of 1 microarcsecond 
(/ias) level. This level of accuracy requires a precise modelling of light propagation. 
In particular, the light deflection due to quadrupole gravitational field of deflecting 
bodies should be taken into account |5J. On the other hand, the accuracy of radio and 
laser radar links of future missions like BepiColombo or Juno require modelling of the 
light travel time at the level of millimeters. Therefore, also the Shapiro delay due to 
quadrupole fields is of practical interest. 

Analytical formulas for quadrupole light deflection are well known. Analytical 
solutions of light deflection in a quadrupole gravitational field have been investigated 
by many authors [H [21 [HI [H [6l [TJ [5l [8] . For the first time the full analytical solution 
for the light trajectory in a quadrupole field has been obtained in [6]. These results 
were confirmed by a different approach in [13]. Various generalization (higher-order 
multipole moments, time-dependence, etc.) were derived in flUi [TT| [T2]. The formulas 
suitable for high-accuracy data reduction are given e.g. in [5]. These formulas are 
rather complicated, so that massive computations of quadrupole light deflection are 
time-consuming. This represents a problem for data processing of astrometric surveys. 
For instance, the ESA mission Gaia (to be launched 2012) will have to process about 
10^^ individual observations of about 10^ distinct celestial objects. It is, therefore, 
obvious that efficient analytical algorithms to compute the quadrupole light deflection 
are mandatory. Furthermore, for observations in solar system, the quadrupole light 
deflection reaches the microarcsecond level only for objects at a relatively small angular 
distance from giant planets. Accordingly, it is highly useful to find simple analytical 
formulas by means of which one can decide whether or not the quadrupole field needs 
to be taken into account for a given accuracy and a given geometrical configuration. 

In this paper we assume that the observer is located within a few million kilometers 
from the Earth's orbit which is clearly the most practical case (e.g. Gaia will have a 
Lissajous-like orbit around Lagrange point L2 of the system Earth-Sun). This allows one 
to simplify the formulas for the quadrupole light defiection considerably. Besides these 
simplified formulas we derive simple analytical estimates of the quadrupole defiection 
allowing one to decide if the effect should be computed and taken into account for a 
given accuracy. We also give a strict upper estimate for the quadrupole Shapiro delay. 

The paper is organized as follows. In Section [2] we summarize some basics about 
light defiection and introduce the notation. In Section [3] the full quadrupole formula in 
post-Newtonian order, the simplest possible expression for stars and quasars still having 
an accuracy of 1 /xas and simple upper estimates of the latter are presented. In Section 
[Hthe full quadrupole formula in post-Newtonian order, the simplest possible expression 
for solar system objects still having an accuracy of 1 fias and criteria are presented. An 
improved estimation of quadrupole Shapiro effect is given in Section [5] The efficiency 
and correctness of the upper estimates and simplified quadrupole formulas have been 
investigated numerically and analytically, and the results are discussed in Section [61 A 
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summary of the findings is given in Section [71 



2. Some basic formulas of light propagation 

Let us summarize some basic formulas of light propagation in post-Newtonian 
approximation. The geodetic equations in post-Newtonian order is linear with respect 
to the metric components and, therefore, the coordinates of a photon and the derivative 
with respect to coordinate time t are given by [5] 

x{t) = x{to) + ca (t - to) + (1) 

i 

x{t) = c(T + Y,^i^i{t) +0{c-^). (2) 

i 

The sum runs over individual terms in the metric of various physical origins (e.g. 

monopole gravitational field of various bodies, quadrupole fields, higher-order multipole 

fields, etc.). Here, to is the moment of emission, Xo = x{to) is the position of source and 
x(t) 

(T = limt^_oo is the unit tangent vector of light path at infinitely past. The position 

c 

of observer is Xi = x{ti) and ti is the moment of observation. The unit coordinate 

direction of the light propagation at the moment of observation reads n = ' | ^\ . In 

\x{ti)\ 

post-Newtonian order the transformation cr to n reads |6l |5] 



n 



<T + ^ 5<Ti + 0{c-*) , dai = (7X (c-^Ai, (ti) X <t) . (3) 



The spherical symmetric part (monopole contribution) due to one massive body A and 
its absolute value are given by (cf. Eq. (102) in [H]) 

KN = -(l + 7)^M.-^(l + ^j. |K.| = (1 + O0j.^(l + ^). 

(4) 

Here 7 is the PPN parameter. Ma is the mass of body A, c is the speed of light, G is 
the gravitational constant, and cIa = cr x (j'f x cr^ is the impact parameter, dA = |<iyi|, 

= x{ti) — xa, = ri , and Xa is the position of massive body A. 
In order to consider light propagation between two given points Xo and Xi (as 
it is needed for the data processing for solar system objects) let us define vectors 
Vq = x(to)—XA and -R = xi—xo = t^ — Tq with absolute values R = \R\ and = |rQ |, 
and unit vector k = R/R. In post-Newtonian approximation, the transformation k to 
n reads [5] 

n = k + J2Ski + 0{c-'^) , Ski = kx(^[c-^Axi{ti)-R-^Axi{ti)]xk). (5) 



The impact parameter cIa can be computed as cLa = k x [rf x kj + 0{c = 

k X (j-Q X k^ + 0{c~'^). The spherical symmetric part (monopole contribution) due 
to one massive body A and its absolute value are given by (cf. Eq. (70) in [5j or 
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Table 1. Numerical parameters of the giant planets taken from [T^. In this table 
the values of are given under assumption that the observer is in the vicinity of 

Earth's orbit. The value J2 for the Sun is taken from [3]. 

Parameter Sun Jupiter Saturn Uranus Neptune 



GMa/c^ [m] 1476. 

[10-3] 0.0002 
Pa [10^ m] 696. 



1.40987 0.42215 0.064473 0.076067 

14.697 16.331 3.516 3.538 

71.492 60.268 25.559 24.764 

0.59 1.20 2.59 4.31 



GM^ Pj/c^ [10^5 m^] 0.143 0.106 0.025 0.000148 0.000165 



Table 2. Maximal numerical values of the neglected terms as given by ([5]) and (|15l) 



Parameter 


Sun 


Jupiter 


Saturn 


Uranus 


Neptune 


[^as] 




< 10- 


^ for Sun and 


giant planets 




[^as] 


1.86 X 10-3 


3.26 X 10 


-2 5.32 X 10" 


-3 8.11 X 10-* 


5.79 X 10-5 



Eq. (24) in [9]) 



-(1 + 7)3- 



rQrf-\-rQ-r 



A ' 



.G Ma 



X ri 



Vq ri + rl 



(6) 



3. The quadrupole light deflection for stars and quasars 

Using the expression Axq (ti) given by Eq. (44) of [5] one gets [HI [17] 



S(T( 



Q ' 



1 + 7 G 
2 c2 



ex 



^1 £a I ^A Va 

A-+^A-+lA— + ^A- 



(7) 



The scalar functions and vectorial coefficients are given by Eqs. flA.ip - flA.SP of 



Appendix A.l The last three terms in Eq. ([7]) can be estimated as 

G Ma 



1 + 7 G 
2 c2 



, , , J^A , Va 

c c c 



p2 
^A 



(8) 



Here, Pa is the equatorial radius, is the second zonal harmonics of massive body A, 
and T-^ IS the minimal distance between massive body and observer. The proof of 
this estimation is given in |17j . 

It is well known (see Table 1 of [5]) that the quadrupole light deflection in solar 
system can achieve the level of 1 /las only for the giant planets (and, possibly, the Sun). 
Using the parameters in Table [1] we obtain from ([8]) that < IQ-^ /ias for all these 
bodies. Furthermore, numerical simulations have conflrmed the correctness of the values 
given in Table [2l Accordingly, these terms in ([8]) can safely be neglected at the level 
of a microarcsecond. Accordingly, the simplest possible expression of quadrupole light 
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deflection for stars and quasars still having an accuracy of 1 /ias and valid for an observer 
situated within a few million kilometers of the Earth orbit reads: 



1 + 7 G 



c 



(9) 



The simphfied formula of quadrupole light deflection (|9]) is still a complicated expression. 
In order to avoid evaluation of this term for each object in the data reduction, a simple 
criterion is needed allowing one to decide whether or not it is necessary to compute the 
quadrupole light deflection for a source. The absolute value of light deflection due to 
the quadrupole field of objects A can be estimated as [17] 



^91 + 7 GMa 




(10) 



A comparison of ( ITOl) with the absolute value of the monopole deflection given by (jl]) 
gives 



9 

< - 



p2 
f_A 

"A 



:iii 



This estimate relates the quadrupole light deflection for stars and quasars to the 
corresponding monopole deflection. The latter is relatively large, defined by a simple 
formula and usually computed for each source and each deflecting body. In this case the 
estimate (fTTj) can be computed at cost of two multiplications (note that d^^ is known 

In case when | 5crpN | is not readily available, one can use 

GMa 



since it is used for 5cr^^ 



<2(l + 7) 



7^ 

'^2 



p2 
f_A 

d\ 



(12) 



<2(l + 7) 



GMa 



Pa 



(13) 



where we use c?^ > Pa and (lA.ip is estimated by |2 + 3 cos a — cos^ a| < 4 for a being 
the angle between vectors cr and r^. Estimate (|T3|) coincides with Eq. (41) of [6]. 



4. The quadrupole light deflection for solar system objects 



The quadrupole light deflection for solar system objects Skq is defined by Eqs. (36)-(47) 
and (69) of [3] and can be written as [HI fTTj : 



6k 



Q 



Y.5k 



A 
Q ' 



7 G 



A 



B 



"a— + Pa — + 7a — + ''a — 



C 



A 

a 



(14) 



A ^ A 

The scalar functions and vectorial coefficients are given in Eqs. flA.lOl) - flA.lTp in 



Appendix A. 2 In [TTj it has been shown that the last three terms in (fT4|) can be 
estimated by 



1 + 7 G 

^2 



< 



2 

'9 



^,1 Ba „ Ca r.,, Va 

Pa — +1a — + ^a — 

c c c 

19 



2 Pir{ 



2 „A min 



+ 



Pa [ri 



A min\ 



+ 



GMa 



^A min\ 
1 ) 



J. 



p2 
^A 



(15) 



On the efficient computation of the quadrupole light deflection 



6 



Using the parameters given in Tabled] we obtain the numerical estimates of given in 
Table [2j Moreover, numerical simulations have confirmed the correctness of the values 
for given in Table [21 In view of these numerical values, the simplest possible form of 
quadrupole light deflection ( !T4l) for solar system objects with an accuracy of 1 fias and 
valid for an observer situated within a few million kilometers of the Earth orbit is given 
by 



6k^ 



1 + 7 G 



A, 



(16) 



2 " c 

Again, in order to avoid unnecessary computations, one needs an efficient way to 
estimate the magnitude of ^fcq. This can be done using the following inequality [T7] 



6k^ 



A comparison with 



< 



3(1 + 7) GMa 1 1 
2 '( 

yields 



"A 



X ri 



A^A^Aj^ ^A 



< 



3n 

2 4 



f)k'^ 



(17) 



This estimate allows one to estimate the magnitude of the quadrupole light deflection in 
observations of solar system objects using the monopole deflection. As it was mentioned 
above, the monopole deflection is significantly larger and should be usually calculated 
for each source and each gravitating body (at least for those bodies, for which the 
quadrupole deflection could be sufficiently large). If | 5fcpN | has been calculated, the 
magnitude of Skq can be estimated at cost of three multiplications (we note that dA is 
required to compute | 6kp^ \ and can be considered as known). If | 6kpj<f \ is not readily 
available, we can use [T7] 

GMa 



<2(l + 7) 



7^ 

■^2 



<2(l+7) 



GMa 



J. 



d\ 
1 

'p'a 



(19) 



(20) 



where in the last estimate we have used Pa < dA- 



5. Shapiro effect for solar system objects 

According to Eq. ([1]), the propagation time cr = c{ti — t^) is given by [6]: 

CT = R + cY,5n + o(c-^) , c5Ti = -k- Axiiti) . (21) 

i 

The formula for the Shapiro delay due to one mass monopole with mass Ma is well 
known: 



C(^T-pN = (1 + 7) — ^ log ^A I ^A 



(22) 

This monopole Shapiro delay becomes unboundedly large for growing distance R 
between the points of emission and observations (although it is growing logarithmically 
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Table 3. Numerical values of estimate (l26l). 



Parameter 


Sun 


Jupiter 


Saturn 


Uranus Neptune 


6 L/2 9 — mm 

1 1 


0.89 


62.16 


20.68 


0.68 0.81 



with R). The quadrupole Shapiro effect c6tq = —k ■ Axqi^ti) is given by [6j: 

1 + 7 G 



(23) 



where the scalar functions and scalar coefficients are given in Eqs. (1A.18I) - (IA.23I) in 



Appendix A.3[ This expression for the quadrupole Shapiro delay cannot be reasonably 
simplified. However, one can give a strict upper bound for the quadrupole effect in the 
Shapiro delay. One can demonstrate that 

GMa 



^\SaVa\<^^ 
c c 

G 

-^\1aJ^a + /3a£a \ < 



p2 
f_A 

^2 ' 



GMa 



Pi 



A\2 



Now since dA> Pa, > Pa, and > Pa, we conclude that 

GMa 



< 3 



•^2 



(24) 
(25) 

(26) 



which represents a strict upper bound of quadrupole Shapiro delay and slightly improves 
the estimate given in Eq. (47) in [6] . This estimate implies that the quadrupole Shapiro 
delay has an upper bound that depends only on physical parameters of the massive 
body. Table [3] gives maximal possible quadrupole effects in the Shapiro delay for any 
positions of the source and observer. 



6. EfHciency of the upper estimates 

As explained above, the principal merit of the simple upper estimates for the quadrupole 
light deflection (Eqs. (fTTj) . (fT2l) and (fT3|l for stars and quasars and Eqs. (fT8|) . (fT9|) and 
(j20l) for solar system objects) is the possibility to use them, at very low computational 
cost, as criteria to decide if the quadrupole deflection should be calculated or not for a 
given configuration and a given numerical accuracy. In this Section we investigate the 
numerical efficiency of the criteria for two situations: (1) purely random homogeneous 
distribution of sources and the position of observer with respect to the deflecting body, 
and (2) light rays grazing the surface of the deflecting body, but with directions still 
randomly distributed with respect to the body (note that the body is not spherically 
symmetric and the orientation does play a role). For both of these situations we compute 
minimal, maximal and mean values of the ratio between the quadrupole deflection and 
its upper estimate. The higher is the mean value the more efficient is the corresponding 
estimate as a criterion. 
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For stars and quasars, starting from fllip . (1121) and (IT^ . we consider the following 
ratios 



ri 



9Pl 
8 d\ 



r2 



•^2 











(a 


•^2 


PI 
d\ 





5(T 


A 
Q 






fA 


7^ 

"^2 


1 
P~A 



(27) 



where ^ctq and 6(7^^^ are determined by Eqs. ([9]) and (jl]), respectively. For solar system 
objects, starting from Eqs. ( fTSl) . ( IT9ll and (j20l) . we consider the ratios 



3Pl 
2 d\ 



r5 





Sk 


A 
Q 






(a 


•^2 







5fc 


A 
Q 








7^ 

•^2 


1 

P'a 



(28) 



where Skq and Sk^^ is determined by Eqs. f|T6|) and ([6]), respectively. 

For all six ratios it is easy to compute the minimal and maximal values analytically. 
Besides that, for both distributions of sources and observers it is possible to compute 
analytically the mathematical expectations (i.e. the mean values) of each of six ratios r^. 
These values are given in Table HI The analytical calculations have been also confirmed 
by direct numerical simulations in which the ratios were computed for correspondingly 
distributed sources and positions of the observer and statistically analyzed. 

The minimal values of all is zero. The maximal values of the ratios are 1 except 
for ri and r4 for grazing rays. In the latter case the maximal values are less than 1. 
This reflects the fact that for grazing rays the numerical coefficients in ( ITTl) and (ITSll 
can be improved (but only for grazing rays and not for arbitrary situation). The fact 
that no maximal values are greater than 1 confirms the validity of the estimates. 

Note that the mean values of and rg for the random distribution depend on the 
ratio Pa/^i of the equatorial radius Pa of the deflecting body and the distance between 
body and observer rf. Although this ratio can be as large as 1, it is small in typical 
applications where the observer is situated far from the body compared to the size of the 
latter. For Gaia the numerical values for the mean values of and rg can be computed 
using Pa and r^™™ from Tabled For example, for Jupiter |(Pa/''i^ ™")^ = 0.49 x 10"^. 

The efficiency of the criteria are characterized by the mean values of the ratios. 
Considering that the random distribution of sources is much more realistic situation 
than the grazing rays we can conclude that the trivial estimates f|T3|) and f l20|) leading 
to r3 and rg are extremely inefficient: the value of quadrupole deflection is typically 
many orders of magnitude lower than "predicted" by those estimates. The criteria (fT2|l 
and ( IT9ll are already better: the quadrupole deflection is typically only 3 times lower 
than "predicted" (the mean value of both r2 and for random sources is 1/3). It is 
clear, however, than the most efficient criteria are given by ( fTTl) and ( fTSl) . For stars 
and quasars the value of the quadrupole deflection "predicted" by ( ITT]) is only two 
times larger than the real value. It means that only 50% of the computations based on 
( ITTl) lead to values lower than the desired numerical cut-off value and could be saved. 
Estimates f lTTj) and f lTSj) will be used for the Gaia data processing. 
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Table 4. Statistical properties of the ratios for two distributions of sources (see 
text for further explanations). 



ratio 




random 






grazing 






min 


mean 


max 


min 


mean 


max 


ri 





40 
81 


1 





16 
27 


8 
9 


r2 
rs 






1 
3 

liPA/rt)' 


' j 


!• 


2 
3 


1 


Ti 





10 
27 


1 





4 
9 


2 
3 


re 






1 
3 

liPA/rty 


' j 


I- 


2 
3 


1 



7. Summary 

In this paper we have developed efficient numerical algorithms allowing one to compute 
the quadrupole light deflection with minimal computational efforts. These algorithms 
will be used for data processing of the ESA astrometric survey mission Gaia and can be 
useful in other cases. In this work we assume that the observer is situated within a few 
million kilometers from the Earth orbit. This is clearly the most practical case. Other 
situations can be analyzed along the lines of our reasoning. The main results which are 
valid with an accuracy of at least 1 yuas are as follows: 

1. Quadrupole light deflection for stars and quasars can be computed as ([9]). 

2. Eqs. ffTT]) . f[T^ and f[T^ can be used as an a priori criterion if the quadrupole light 
deflection ([9]) has to be computed for a given source. 

3. Quadrupole light deflection for solar system sources can be computed as f|T6|) . 

4. Eqs. ( |T8l) . ( |T9l) and ( l20l) can be used as an a priori criterion if the quadrupole light 
deflection (fT6l) has to be computed for a given solar system object. 

The efficiency of the upper estimates has been investigated numerically and analytically, 
and the results are shown in Table |H They demonstrate high efficiency and correctness of 
the upper estimates, both for randomly distributed sources and sources which generate 
grazing rays. According to these investigation, the most efficient upper estimate of 
quadrupole light deflection is f lTTj) for stars and quasars and f|T8|) for solar system objects. 
The correctness of simplified quadrupole formulas has also been shown by numerical 
simulations. 

Additionally, we give a strict upper bound (!26l) for the quadrupole effect in the 
Shapiro delay. This upper bound can be used to decide whether or not the quadrupole 
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Shapiro delay should be taken into account if high-accuracy ranging measurements are 
to be modelled; e.g. in the framework of missions like BepiColombo [15] or Juno. 
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Appendix A. Explicit formulas for the quadrupole terms 

We give the full expressions of coefficients and scalar functions of quadrupole light 
deflection and quadrupole Shapiro effect in post-Newtonian order, because so far they 
were not being presented in a refereed journal. 

Appendix A.l. Light deflection for stars and quasars 
The functions in ([7]) read 

^-:^f2 + 3^-i^|, (A.1) 
c d\y rt {r^f ' ' 



(A.2) 



(A.4) 



c {rt) 
Va ^ \_ 

a^' = -M^ g + 2 M^. ^ - 2 aV^ ^ - 4 , (A.5) 

/5? = 2M>^^, (A.6) 

Ji jj Jk Jk 
J k l\/rA "A "A "A ,i^A i 1- "A ( \ ^\ 

^ — 4 Ta' ^ ' 

6'/ = -2 M,j a' a'' + 2 M^^ - AMf^a' . (A.8) 

"A 

Here Mf- is the symmetric and trace-free quadrupole moment of body A. For an axial 
symmetric body (this approximation is sufficient for the solar system and the accuracy 
of 1 /ias) one has 

Mf^ = ^MA Pi 4 T^A diag (1, 1, -2) , (A.9) 

where TZa is the rotational matrix giving the orientation of the figure axis of body A 
(see Eq. (48)-(53) of [5J). 
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Appendix A. 2. Light deflection for solar system objects 
The functions in ffT4l) read 



c dA R \rQ r^ — k ■ Tq r^r^ — k-rij d\ \ r^ 



Ba _l (k-r^ k-r^\ ^ [rf) - 3 (fc ■ rf 



c 



R\{r^f {r^fl {rtf 

CA_dA(j^ l_\ _ k-rf 

c R[ir^f (rtf) ''\rtf' 
Va 1 1 fk-r^ k-rf\ 1 



c 



d\ R \ rt ) {rtf ' 



^ dA ^ dA ^ dA ^ d\ 

P'^' = 2M:^ k^ ^ , 

Ua 

ji jj jk jfc 
" ^^'^ d\ ^^'^ ^ ^ dA' 

5'^^ = -2Mf^k'k^k'' + 2M^^k^ - AM^jk'^^. 

a /i 



Appendix A. 3. Shapiro delay 

The functions in (123|) read 
k ■ k ■ 



1 ( k ■ k ■ rf' 



Va = -^ 



d\ \ rff ri J 

r/* d^ 

"A "A 
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